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Abstract. Closed-loop or feedback controlled ratchets are Brownian motors that operate using information 
about the state of the system. For these ratchets, we compute the power output and we investigate its 
relation with the information used in the feedback control. We get analytical expressions for one-particle 
and few-particle flashing ratchets, and we find that the maximum power output has an upper bound 
proportional to the information. In addition, we show that the increase of the power output that results 
from changing the optimal open-loop ratchet to a closed-loop ratchet also has an upper bound that is 
linear in the information. 
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1 Introduction 

Brownian ratchets have been studied in different contexts 
due to their theoretical importance in non-equilibrium sta- 
tistical mechanics and their potential relevance for appli- 
cations in disciplines like nanotechnology, condensed mat- 
ter or biology [HdJEll!] • Many studies deal with the per- 
formance of these devices (see refs. [TJ[5] for comprehen- 
sive reviews) concentrating on open- loop ratchets, as those 
obtained fluctuating an uniform external force (rocking 
ratchets [6,7 J, or an external asymmetric potential (flash- 
ing ratchets [6,8J), either randomly or periodically. On 
the other hand, closed-loop or feedback controlled ratch- 
ets, as the so-called instant maximization protocol [9] and 
the threshold protocol [TUj, use information of the state 
of the system to operate. The feedback ratchet of [9] has 
been recently proposed as an effective model to describe 
the stepping motion of the two-headed kinesin [3] . Other 
'information-dependent' rectification mechanism have been 
recently proposed to model certain chemical and biological 
systems [4]. 

The previous works [9l [T0lfTT] about closed-loop ratch- 
ets focussed on the study of the flux and its maximization. 
In particular, it has been shown that the increase of the 
flux performance when the optimal open-loop control is 
changed to a closed-loop control has an upper bound pro- 
portional to the square root of the information used by 
the controller ITTI . 



In this paper, we consider another measure of the per- 
formance, viz. the power output, with the aim of getting 
further insight in the relation between information and the 
increase of performance in a system with thermal fluctu- 
ations. We oppose to the flux a constant load force [12] in 
order to compute the potential energy gain by the parti- 
cles thanks to the action of the controller. The generaliza- 
tion of the methods developed in [TT] allow us to obtain 
the relations between the maximum power output and the 
information that the controller uses. 



2 The model 

The collective feedback ratchet that we investigate has 
two basic ingredients, namely, N Brownian particles and 
a controller. The controller acts on the particles switch- 
ing on and off a potential V(x) according to the control 
policy and to the information received about the state of 
Brownian particles through a noisy channel. 

Specifically, we consider N overdamped Brownian par- 
ticles at temperature T in a piecewise linear saw-tooth 
potential 



V(x) = 



&(f-*) 



if < f < a, 
if a < f < 1, 



(1) 
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of height Vq, asymmetry parameter a, and period L, i.e. 
V(x) = V(x + L). The potential is switched on and off 
according to the instant maximization of the center-of- 
mass velocity protocol (see ref. [9]), which switches on 
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the potential only when the net force due to the poten- 
tial on the particles would be positive. In order to obtain 
work from the system operation we oppose to the flow 
of particles an homogeneous static force F ext ', thus, the 
total force acting on the particles when the potential is 
on is F tot (x) = F(x) - F ext , with F(x) = -V'(x), and 
Ftot(x) = —F cxt when the potential is off. The state of 
the system is described by the positions Xi(t) of the par- 
ticles that satisfy the Langevin equations 

7x»(t) = a(t)F( Xi (t)) - + &(*); i = l,...,N, (2) 

where 7 is the friction coefficient (related to the diffusion 
coefficient D through Einstein's relation D = kuF '/j) and 
£i (t) are Gaussian white noises of zero mean and variance 
(O) = 1lk B T8ij8{t - t'). The dichotomous func- 
tion a(t) [a = (potential off) or a = 1 (potential on)] 
implements the action of the controller. The control pol- 
icy uses the information received from the system through 
a noisy channel that we model with a binary symmetric 
channel [13] . This channel passes the sign of the net force 



1 N 



(3) 



to the controller with an error probability p known as the 
noise level of the channel, so when f(t) > (< 0) the 
controller switches on (off) the potential with probabil- 
ity 1 — p. Therefore, the feedback protocol and the noisy 
channel lead to the effective control policy 



a eS (t) = (l-p)G(f)+pO(-f), 



(4) 



with the Heaviside function [0(x) — 1 if x > 0, else 
0(x) — 0]. This effective control policy is equivalent to the 
protocol of instant maximization through a noisy channel 
provided many measurement and control actions are per- 
formed in the characteristic time of the system evolution, 
which is the case we consider here. 

Our aim is to study the dependence of the maximum 
power output with the information. On one hand, the av- 
erage information transmitted through the noisy channel 
is quantified in terms of the mutual information 13J that 
the controller gets from the state of the system. Our case 
— the noisy measurement of the sign of the net force — is 
equivalent to a noisy channel called the binary symmetric 
channel in information theory. For this case the mutual 
information can be computed (see sec. 8.1.4 of [13]), and 
it is given (in bits) by 

I = H(q)-H(p), (5) 

with H(x) = ~x\og 2 x — (1 — a;)log 2 (l — x), q = (1 — 
p)b + p(l — b) the probability that the controller receives 
a negative sign, and b the probability that the actual sign 
of the net force is negative. Therefore, the information I 
that the controller gets about the system is greatly deter- 
mined by the noise level p of the channel; the maximum 
information is reached for p — and it is at most 1 bit, 
while for p — 1/2 the channel becomes completely ran- 
dom and no information of the system is received by the 



controller. When the probability b does not depend on p, 
the relation between the noise level p and the information 
/ can be easily expanded around p = 1/2 and reads 



'M = i^£^[i-<i-»)'](p-5 



(6) 

Inverting this relation to leading order we get for p < 1 /2 
the result [TT] 



/In 2 



(1-6) 



(7) 



On the other hand, a positive power output is obtained 
when there is a net flux against the load F cxt that tilts 
the potential. In the stationary regime the center-of-mass 
moves with a mean velocity (x cm ) and then the average 
power output (work obtained per unit time) is given by 



F F cx t(x cm ) 



(8) 



We first analyze the dependence of the power output 
with the information for the case of one particle and later 
for the few-particle ratchet. 



3 One particle 



We start with the one particle case (N = 1) where an 
effective potential that includes the effects of the load, 
the control protocol, and the binary symmetric channel 
can be constructed. The system dynamics can be viewed 
as the result of the action of the periodic effective force 
F e g(x) — a c g(x)F(x) — F CKt that derives from an effective 
potential. Using units L = 1 and fc^T = 1, this effective 
potential can be written in terms of K :— pVo + F cxt a and 
M:= (1 - p)V Q - F cxt (l - a) as 



V cS (x) = 



xK 

K — Mf^ 2 

l-a 



if < x < a, 
if a < x < 1 



(9) 



in the interval [0, 1], and outside V c s(x) = V c g(y) + (x — 
y)Ve ff (l), with y = x mod 1, y € [0, 1]. eq. © and fig. [T] 
show that the effect of increasing the noise level p is to 
diminish in the effective potential the average tilt that 
induces a positive flux, while the effect of increasing the 
load F cxt is to tilt the effective potential opposing the 
positive flux. 

Solving the stationary Fokker-Planck equation for this 
effective potential the stationary mean velocity for one 
particle is obtained (in units L — 1 and ksT = 1): 



(x) = 



DK 2 M 2 A 
AE-B+B- 



(10) 
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Fig. 1. Effective potentials for one particle and noise levels p = 
0, 1/4, and 1/2 with Vo = 5 and a = 1/3. Left panel without 
external load, and right panel with external load F cxt = 1. 
Units: L = 1, D = 1, and k B T = 1. Note that for p = 1/2 the 
effective potential is equal to V(x)/2 + F cxt x. 




Fig. 2. Power output versus the load for Vo = 5 and a = 1/3 in 
the one particle case [eqs. ((8} and (|10|l ]. Units: L = 1, D = 1, 
fc B T = 1. 



with 



A := 
B± := 

J5 := 



[aM + (1 - a)if]e ±K 
- (1 - a)Ke^ K - M ^ - aM, 
a 2 M 2 (l -K -e~ K ) 



+ a(l - a)KM(l - e M )(l 
+ (1 - a) 2 K 2 (l + M ■ 



For p < 1/2, there is a positive current for forces 
smaller than the "stopping force" F st0 p (the one that leads 
to the cancellation of the velocity), so a work is done 
against the load for F ext € (0, F s t op ). For p — 1/2, the 
stopping force is zero, because no positive flux is obtained 
even in the absence of the external load. Our noisy con- 
trol acts instantaneously, i.e. in a time scale much faster 
than the characteristic times of the system \(aL) 2 /{2D) 
for the diffusion time and 7(1 — a) 2 L 2 /Vq for the charac- 
teristic time of the drift induced by the potential]. Thus, 
for p = 1/2 the potential V(x) is randomly switched on 
and off very fast and the particle just feels the average po- 
tential. This implies that the effective potential in absence 
of the load, namely V(x)/2 (see fig. [T|) , is not tilted, giving 
a zero flux for the p = 1/2 case for zero load. Therefore, 
in order to get work the noise level of the channel should 
be p € [0, 1/2) and the value of F stop is obtained equating 
eq. (jTUJ) to zero, 



^ t op = ^(l-2p). 



(11) 



Substituting eq. (fTI))) in eq. ([5]) we get the analytical ex- 
pression for the power output in the one-particle ratchet. 
The dependence with the load is plotted in fig. [2] for noise 
levels p = 0, p = 1/4, and p = 1/2. The positive regions 
correspond to the system doing work against the external 
force. The load F* xt that maximizes the power output lies 
between and F stop and it is given by the condition 

dP 



dF c: 



0. 



(12) 



In general, it is a function of the noise level of the binary 
symmetric channel and it also depends on the physical 
parameters of the potential, Vq and a. The condition (fT"2l) 



gives a transcendental equation for F* xt that can be nu- 
merically solved in order to obtain the maximum power 
output, 



(13) 



This equation gives the dependence of the maximum power 
with the noise level p, which is related with the informa- 
tion / through eq. ([5]) . This last equation requires to com- 
pute b, which can be obtained integrating over the space 
interval [0, aL] the stationary distribution of the Fokker- 
Planck equation for the effective potential ([9]), 




-K\ 



1 - e- K + (1 - <z)(l 



-M 



)K/{aM) 



-K 



-M 



(14) 



-K 



(units L = 1, D = 1, ksT = 1). Therefore, the combi- 
nation of eqs. ©, and (|12H14p permits to obtain the (im- 
plicit) exact dependence of the maximum power developed 
by the Brownian motor as a function of the information 
gathered by the controller (see fig.[3|). 

We analyze now the regime of small potentials in the 
one particle case. For small potentials (Vo < ksT) the 
value of the external force that maximizes the power is 
also small [remember eq. (fTTj) and the fact that F* xt < 



(15) 



-Fstop]- In this regime, the velocity (fT0|) reduces to 
D(M — K), or, recovering units, 



(x) ~ (1 - 2p) 



7 



We see that there are two contributions to the velocity: 
the current effect due to the white thermal noise and the 
control through the binary channel, (1 — 2p)Vo/ (^L) : and 
the net drift due to the load, —F ext /^- We highlight that 
for small potentials (and loads) these two effects appear 
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Fig. 3. Maximum power output in the one particle case as a 
function of the information used. The curves are for heights of 
the potential Vq = 2, 5, 10 and asymmetry parameter a — 1/3. 
Units: L = 1, D = I, k B T = 1. 
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Fig. 4. Maximum power output in the one particle case as 
a function of the information for Vo = 1 and a = 1/3, and 
comparison with the upper bound (|18[1 and the better upper 
bound ([20j). Units: L = 1, D = 1, k B T = 1. 



uncoupled, and the result is independent of the asym- 
metry of the potential. This independency of the asym- 
metry for small potentials can be understood realizing 
that in this case the effective potential is well approxi- 
mated by a flat potential with the same average slope, 
i.e., V eff {x) ~ [-Vb(l - 2p)/L + F cxt ]x. 

Applying eq. (|12[) to the power output computed using 
cq. (fT5T) we obtain 



2L 



(1 - 2p) 



F, 



stop 



and then the power output is 



F *2 



V 2 

(l-2p) 2 . 



4 7 L 2 



(16) 



(17) 



On the other hand, for small potential heights b ~ a, and 
using eq. ([7]) we get 



Ril, 



(18) 



with Ri a constant that depends on the physical param- 
eters of the system, 



Ri 



Vo 2 ln2 



87L 2 a(l -a)' 



(19) 



Notice that the dependence on the asymmetry a does ap- 
pear here because it determines the relation between p 
and I [eq. (0)], as b ~ a for small potentials. 

Therefore, eq. (fT7|l indicates that for small potential 
heights and small values of the information (i.e., p ~ 1/2) 
the maximum power is approximately directly propor- 
tional to the information gathered. In addition, we have 
numerically checked that eq. (|18p gives an upper bound of 
the maximum power for any potential height Vq and for 
any value of the information /. 

A better approximation for the dependence of the max- 
imum power output with the information can be found 
using the result of inverting eq. ([6]) up to fourth order, 



P„ 



Si(-1+ y/l + S 2 I), 



(20) 



with 



Si = 



w 2 



l-(l-2a) 2 



5 2 = - 



4 l-(l-2a) 4 



■In 2, 



4 7 L 2 1 - (1 - 2a) 4 ' " 3 [1 - (1 - 2a) 2 ] 2 

(21) 

which is also an upper bound of P max for any potential and 
information values. In fig.Uthese upper bounds [eqs. (|18p 
and (f2"0")) ] are compared with the exact result for Vo = 
k B T. 



4 Few particles 

Let us now study a collective ratchet composed of a few 
particles and show that the results are similar to those in 
the one particle case. Summing and averaging the Langevin 
equations (|2|), the average velocity of the center-of-mass 
in the stationary state can be written as 



7<icm) = (a e ff(/)/) - F c: 



(22) 



An approximate solution can be found assuming, as in |11[ 
|9], that: (i) the position of the particles are statistically 
independent, and (ii) the probability of finding a parti- 
cle in a negative force interval (for example [0, aL]) is a. 
These assumptions are verified for small potentials and 
small loads even in the presence of noise, and they imply 
that the probability distribution for / is approximately 
Gaussian, 

with £ the amplitude of the fluctuations of the net force, 
given by S = — . Vo (see ref. fSl). Following [TT], we 

LWa(l—a)N 

get for the center-of-mass velocity 



^ , v Ff,vf 

— =(1 - 2p) — . 

7V2T 7 



(24) 



This result is the sum of the center-of-mass velocity with- 
out the external load [TTj plus the drift —F cxt /j due to the 
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external load. We see that, like in the one particle case, 
these two effects are decoupled for small potential heights. 
The expression (|24l) agrees with the results of numerical 
simulations of the stochastic evolution equations {2} . 

Applying eqs. ([8]) and ([12")) . it can be shown that the 
maximum power is reached for 



V 



Lyj8Tra(l-a)N 



(1 - 2p) 



and takes the value 



P, 



V 2 

7L 2 8vra(l - a)N 



(l-2p? 



(25) 



(26) 



or simply P max = F* 2 t /j. Therefore, we also have in the 
regime of few particles (and small potentials) that the 
maximum power is at first approximation directly pro- 
portional to the information, 



P, 



Rjyl 



(27) 



In the previous expression, the constant Rjy depends only 
on the physical parameters of the system, in particular the 
number of particles N, 



Rn = 



Vo 2 ln2 



(28) 



7 L 2 167ra(l - a)6(l - b)N' 
where b can be calculated for small potentials and loads 



JV 
n>aN 



N 



a n (l-a) 



N-r 



(29) 



using the same assumptions that lead to the Gaussian ap- 
proximation for /?(/). We have checked numerically that 
eq. (I27p is an upper bound for the maximum power out- 
put in the few particles case. Again, as in the one-particle 
ratchet, a linear upper bound has been found for the maxi- 
mum power output that the system can get using a certain 
amount of information I. 



5 Correlation 

The previous expression of the different physical magni- 
tudes of the system in terms of the noise level p can be 
recast in terms of a correlation C that we introduce in this 
section. The main underlying idea is that the presence 
of noise in the control induces a decorrelation between 
the relevant magnitudes of the control policy. In the in- 
stant maximization of the center-of-mass protocol [9] the 
switching of the potential only depends on the sign of the 
net force, namely sgn/, and the presence of noise in the 
control implies that the controller does not use the ac- 
tual value sgn/ but a value sgn/. The correlation between 
these quantities, 



C = (sgn/ sgn/), 



(30) 



can be written as 



C = P + 



P__ - P, - P- 



(31) 



where P^ is the join probability of having sgn/ = +1 

while the controller receives sgn/ = — 1, and analogously 
for the other joint probabilities. As in our system sgn/ 
is different from sgn/ with probability p (the noise level) 
these joint probabilities can be easily computed by noting 



that P_ 



+ 



bp, P- 



6(1 - p), P + 



(1 — b)p and 



P ++ = (1 — b)(l — p), with 6 the probability of sgn/ being 
— 1. Therefore, the correlation can be parameterized in 
terms of the noise level p as 



C=l-2p. 



(32) 



In other words, the effect of the noise is to decrease the 
correlation C, which has its maximum value (C = 1) for 
zero noise and its minimum (C = 0) for a completely noisy 
policy, p = 1/2. 

Finally, using eq. (|32|) . the relations derived in previous 
sections can be restated in terms of the correlation. For 
example, eqs. (fTT)) and (|2"6"1) reads 



(one particle case), and 
p = 

1 max 



V a 2 C 2 
4-fL 2 ' 



V 2 C 2 



7L 2 87ra(l-a)A^' 



(33) 



(34) 



(few particles case). 

This reformulation helps to understand the physical 
meaning of the relations derived in the previous sections 
giving a complementary view. In addition, it indicates that 
the noisy control considered can give an effective descrip- 
tion of other feedback ratchets with an imperfect opera- 
tion of the feedback control. For instance, this effective 
description has been shown to be valid in time-delayed 
feedback ratchets consisting of few particles [H], where 
/ = fit), f = f(t — r) (with r being the time-delay), and 
the correlation can be computed just from the time series 
of the net force f(t). 



6 Comparison with open-loop protocols 

The instantaneous maximization of the center-of-mass ve- 
locity is the optimal protocol to maximize the power out- 
put in the one particle case for a noiseless channel (p — 0). 
Thus, we also expect this protocol to give a power output 
close to the maximum possible value in the few particles 
case and in the presence of noise with a memoryless pro- 
tocol (note that protocols with memory can perform error 
corrections) . Therefore, we expect eq. (IT8|) and eq. (|27|) to 
be upper bounds of the maximum power output that can 
be obtained with a memoryless closed-loop protocol that 
uses an amount of information / about the system. 

In addition, the maximum power output obtained with 
open-loop protocols is much smaller than that obtained 
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with efficient closed-loop protocols. For instance, for the 
saw-tooth potential with parameters Vq = 5&bT and a — 
1/3, the periodic protocol with optimum periods T on ~ 
0.06L 2 /D and T g — 0.04L 2 /D gives a small maximum 
power -P^p™ — 0-04V r 2 /(7L 2 ), which is reached for a load 
F* xt ~ 0.25Vo/i- In contrast, the closed-loop one-particle 
ratchet yields a maximum power P^^ cd — 5.1Vq / ('yL 2 ) 
for F* xt ~ 2AVq/L when it works without noise in the 
channel. Therefore, the linear eqs. (fT8|) and (|27|) are also 
good estimates of the maximum improvement that can be 
attained changing from the optimal open-loop control to 
a closed- loop protocol, i.e., 

pclosod _ popon / nr /or\ 
1 max 1 max — Ju 3 \ 00 J 

where R is a constant depending on the system's charac- 
teristics; see eqs. (fl9|) and ([28]). 



7 Concluding remarks 

In this article we have analyzed the relation between the 
information about the state of the system used by the 
controller and the power output in a feedback controlled 
ratchet. We have obtained exact analytic results for one- 
particle ratchets, and also approximate simple expressions 
for the maximum power output in both one-particle and 
few-particle ratchets. Moreover, we have found that the 
increase of the maximum power output when we change 
from the optimal open-loop protocol to a closed-loop pro- 
tocol has an upper bound proportional to the information 
used by the controller. Also an upper bound proportional 
to the information was found in [15j for the entropy reduc- 
tion in a general closed-loop controlled system. The result 
obtained in the present paper for the maximum power out- 
put is the analog upper bound of the one found in [11| for 
the flux, but with the important difference that the upper 
bound for the flux was proportional to the square root of 
the information. 
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